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A  Composite  Thermoelectric  Device  (CTED)  is  comprised  of  thermoelectric  (TE)  elements  made  of  TE 
semiconductor  materials  bonded  to  highly  electrically  and  thermally  conductive  material  in  a  segmented 
fashion.  Thermoelectric  performance  of  such  devices  using  numerical  methods  with  temperature- 
dependent  thermoelectrical  properties  has  been  investigated.  The  CTED  performance  in  terms  of 
produced  electrical  current,  Ohmic  and  Seebeck  potentials,  power  output  P0,  heat  input  Qh,  and 
conversion  efficiency  77  is  studied  for  various  hot  surface  temperature  Th,  load  resistance,  semiconductor 
thickness  d,  and  convection  heat  transfer  coefficient  h  values.  The  aforementioned  CTED  performance 
characteristics  are  compared  to  those  of  a  conventional  TED  with  geometrical  equivalence.  For  a  given  Th, 
a  maximum  Po  is  achieved  at  a  load  resistance  value  that  is  equal  to  the  total  internal  resistance  R[  of  the 
device;  an  optimum  77,  at  an  optimum  load  resistance  R0ptmL  which  is  typically  higher  than  the  Ri.  A  CTED 
with  d  =  1  mm,  the  optimum  77  values  are  24.8%,  26.2%  and  29.9%  higher  than  conventional  TED  values  at 
Th  =  350  K,  450  K  and  550  K,  respectively.  At  R0ptmL  values,  the  difference  in  P0  and  Qh  show  significant 
and  minor  increases,  respectively,  in  relation  to  differences  in  77  with  an  increase  in  Th.  The  variation  of 
the  semiconductor  thickness  d  has  a  substantial  effect  on  the  CTED  characteristics  and  R0ptmL  values;  as 
d  decreases,  a  continuous  increase  in  Po  and  Qh  and  an  optimum  value  of  77  are  achieved.  Intuitively, 
PoptmL  increases  with  an  increase  in  d  and  it  reaches  a  maximum  value  at  the  conventional  TED  limit. 
With  d  =  0.5  mm,  Th  =  450  K  and  h  =  20  W  m-2  K-1  at  a  corresponding  Roptmi  value,  P0  and  Qh  exhibit  an 
eight-  and  six-fold  increase,  respectively,  and  77  is  increased  22%  compared  to  a  conventional  TED.  The 
convective  heat  transfer  coefficient  has  a  pronounced  effect  on  the  CTED  performance  when  it  is  greater 
than  100  W  m-2  K_1.  From  this  study,  the  CTEDs  show  promise  of  extracting  more  heat  in  waste  heat 
recovery  applications  when  compared  to  conventional  TEDs. 

©  2012  Elsevier  Masson  SAS.  All  rights  reserved. 


1.  Introduction 

Approximately  two-thirds  of  the  fossil-fuel  energy  we  use  in 
automobiles,  power  plants,  and  industrial  systems  is  rejected  into 
the  atmosphere  as  waste  heat.  A  viable  technology  for  efficiently 
converting  waste  heat  into  electricity  is  thermoelectric  devices 
(TEDs).  Based  on  the  principle  of  Seebeck  effect,  these  devices  work 
as  a  power  generators  by  maintaining  a  temperature  difference 
between  hot  and  cold  junctions.  They  have  no  moving  parts,  emit 
no  noise  and  are  environmental  friendly. 

Over  the  past  few  decades,  TEDs  made  of  bulk  bismuth- 
antimony—  telluride— selenium  (Bi,  Sb,  Te,  Se)  have  not  shown  an 
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increase  in  either  their  performance  (Figures  of  Merit  (ZT)  ~  1)  or 
thermal  conversion  efficiencies  (~5— 10%)  [1,2].  However,  with  the 
growing  environmental  consciousness,  there  is  a  renewed  interest 
in  enhancing  these  renewable  and  sustainable  device’s  perfor¬ 
mance  for  electrical  power  generation  and  waste  heat  recovery 
applications  (up  to  ~30%)  via  methods  of  nano-structuring,  nano¬ 
fabrication  and  use  of  new  bulk  materials  [3].  Caillat  et  al.  [4] 
developed  a  segmented  TED  which  operates  over  a  temperature 
range  of  300-973  K  using  novel  p-type  (Zn4Sb3,CeFe4Sbi4)  and  n- 
type  (CoSb3)  based  alloys.  They  achieved  a  conversion  efficiency  of 
15%  and  also  demonstrated  the  bonding  studies  between  the 
segments.  Gou  et  al.  [5]  investigated  the  performance  of  low 
temperature  waste  heat  recovery  TEDs  using  one-dimensional  (ID) 
numerical  simulations  and  experiments.  They  concluded  that  in 
addition  to  increasing  the  waste  heat  temperature  and  number  of 
modules  in  series,  expanding  the  heat  sink  surface  area  and 
enhancing  the  cold  side  heat  transfer  in  proper  ranges  can  have 
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Nomenclature 

£ 

emissivity,  dimensionless 

V 

thermal  efficiency,  dimensionless 

A 

cross-sectional  area,  mm2 

K 

thermal  conductivity,  W  m-1  I<-1 

As 

surface  area,  mm2 

P 

electrical  resistivity,  Q  m 

d 

semiconductor  thickness,  mm 

O 

Boltzmann  constant,  W  m2  I<4 

D 

depth,  mm 

h 

convection  heat  transfer  coefficient,  W  m-2  K-1 

Subscripts 

H 

height,  mm 

c 

cold  wall,  conductor 

I 

electric  current,  A 

i 

internal 

J 

electric  current  density,  A  m“2 

ic 

inter-connector 

J 

electric  current  density,  vector 

h 

hot  wall 

L 

length,  mm 

L 

load 

Po 

power  output  (J2RL),  W 

n 

n-type  semiconductor 

q" 

heat  flux,  W  m“2 

optmL 

optimum  load 

a 

heat  transfer,  W 

0 

Ohmic 

T 

temperature,  K 

P 

p-type  semiconductor 

V 

voltage,  V 

S 

Seebeck 

R 

electric  resistance,  Q 

sur 

surroundings 

W 

width,  mm 

? 

normal  surface 

00 

ambient 

Greek  symbols 

a 

Seebeck  coefficient  or  thermopower,  V  I<  1 

dramatic  effects  on  TED’s  performance.  Furthermore,  Hsiao  et  al.  [6] 
also  developed  a  1 D  thermal  resistance  model  for  thermoelectric 
modules  with  applications  to  waste  heat  recovery  from  an  auto¬ 
mobile  engine  and  compared  their  model  results  with  experi¬ 
mental  data.  Their  outcomes  showed  that  the  thermoelectric 
module  on  the  exhaust  pipe  performs  better  compared  to  a  module 
on  the  radiator  system. 

Most  of  the  studies  in  literature  were  investigated  with  ID 
models,  even  though  electric  field  and  temperature  distributions 
are  multidimensional  in  TEDs.  Moreover,  the  heat  flux  and  current 
density  are  not  in  parallel  in  conventional  TEDs.  Therefore, 
considering  the  three-dimensional  (3D)  effects  in  numerical 
models  will  give  more  realistic  characteristics  of  TED  performance. 
It  has  been  shown  in  the  articles  [7-15]  that  3D  simulations  are 
more  accurate  in  design  and  optimization  of  TEDs.  They  conducted 
numerical  simulations  of  TEDs  with  in-house  codes  as  well  as  with 
available  commercial  software  packages.  Harris  et  al.  [7]  developed 
a  finite  volume  based  numerical  code  with  temperature-dependent 
material  properties  to  study  3D  temperature  and  electrical  field 
distributions  in  TEDs.  They  investigated  the  device  performance  by 
taking  into  consideration  both  thermal  and  electrical  contact 
resistances  at  the  interfaces  of  the  thermoelectric  material  and 
conducting  shoes.  Additionally,  they  analyzed  the  TED  with  inert 
gas  and  insulating  materials  between  element  legs,  to  reduce  and 
eliminate  convection  between  hot  and  cold  surfaces  within  said 
device,  respectively.  Chen  et  al.  [8]  analyzed  the  thermoelectric 
generator  performance  connected  to  various  models  of  load 
circuitry  in  a  SPICE-package  environment  and  validated  their 
results  with  ANSYS  and  experimental  results.  Further,  Hu  [9] 
studied  the  characteristics  of  two  compact  gas-phase  heat 
exchangers  placed  on  the  hot  and  cold  sides  of  the  TED  to  enhance 
heat  transfer  and  fluid  flow  while  minimizing  thermal  stresses. 

Ziolkowski  et  al.  [10]  analyzed  the  performance  of  thermo¬ 
electric  generators  including  the  parasitic  effects  such  as  convec¬ 
tion,  radiation  and  conduction  heat  bypass.  They  also  considered 
the  design  parameters  such  as  pellet  aspect  ratios  and  thermal  and 
electrical  contact  resistances.  Kousksou  et  al.  [11]  examined  the 
thermoelectric  power  generation  from  a  helicopter  conical  nozzle. 


They  showed  that  the  electric  power  produced  is  sensitive  to  the 
hot  fluid  properties  and  significantly  depends  on  actual  operating 
conditions.  Using  Technology  Computer-Aided  Design  (TCAD) 
package,  Gould  et  al.  [12]  simulated  the  thermoelectric  couple 
configured  for  low  power  generation.  Chen  et  al.  [13]  proposed  and 
implemented  a  3D  numerical  model  for  thermoelectric  generators 
in  FLUENT  UDS  (User  Defined  Scalar)  environment.  This  model 
accounts  for  all  temperature-dependent  properties  of  materials  and 
non-linear  fluid-thermal-electric  multiphysics  coupled  effects.  The 
model  was  validated  with  simulation  data  from  other  models  and 
experimental  data.  Recently,  in  the  articles  [14-16],  using  3D 
numerical  methods  implemented  in  FLUENT  UDS  package  and 
analytical  solutions,  the  authors  studied  the  enhanced  heat  transfer 
characteristics  and  performance  of  composite  and  integrated  TEDs 
applied  to  waste  heat  recovery  systems.  Further,  various 
researchers  studied  optimum  geometric  dimensions  to  achieve 
maximum  conversion  efficiency  of  TEDs  [17,18]. 

A  review  of  literature  suggests  there  is  increasing  interest  in 
advanced  TE  materials  research  for  higher  ZT  and  the  exploration  of 
novel  TED  designs  for  increased  performance  while  promoting  the 
notion  of  sustainability.  State-of-the-art  and  sustainable  TED 
designs  constructed  with  minimum  amounts  of  rare-earth  material 
without  compromising  conversion  efficiency  and  performance 
while  increasing  the  power  output  is  one  of  the  methods  to  mini¬ 
mize  the  device’s  cost.  A  potential  and  promising  approach  is  to  use 
a  novel  composite  thermoelectric  device  (CTED),  where  a  major 
part  of  the  semiconductor  is  replaced  with  a  conductor  such  as 
copper.  In  this  work,  by  concentrating  on  the  structuring  of  bulk 
materials,  the  performance  of  CTEDs  compared  with  conventional 
TEDs  has  been  investigated  using  a  3D  numerical  model  imple¬ 
mented  in  the  FLUENT  UDS  environment. 

2.  Geometry,  mathematical  model  and  boundary  conditions 

The  three-dimensional  schematics  of  composite  and  conven¬ 
tional  thermoelectric  devices  (TEDs)  being  investigated  are  shown 
in  Fig.  la  and  b,  respectively.  The  TED  consists  of  two  vertical  legs 
connected  electrically  in  series  and  thermally  in  parallel  via 
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Composite  thermoelectric  device  Conventional  thermoelectric  device 
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Fig.  1.  Schematic  of  (a)  composite  and  (b)  convectional  thermoelectric  devices. 


connectors  made  of  a  highly  electrical  conducting  material.  The 
composite  TED  (CTED)  leg  consists  of  a  copper  inter-connector 
interposed  between  bulk  crystalline  n-type  or  p-type  materials, 
whereas  a  convectional  TED  leg  is  entirely  made  of  either  n-type  or 
p-type  material.  The  semiconductor  thickness  d  has  been  taken 
equal  in  size  for  both  the  left  and  right  legs  of  the  CTED.  Each  leg  has 
a  square  cross-section  WxD  and  height  H,  and  the  legs  are  separated 
by  a  distance  L\  as  shown  in  Fig.  1.  The  top  and  bottom  surface  of  each 
leg  are  connected  to  copper  connectors.  The  end  surfaces  of  these 
connectors  act  as  electric  terminals  for  load  circuitry  of  load  resis¬ 
tance  (Kl).  The  top  surface  of  the  upper  connector  is  subjected  to 
constant  hot  temperature  (Th)  while  the  bottom  surface  of  the  lower 
connector  is  maintained  at  a  constant  cold  temperature  (Tc).  The 
remaining  surfaces  of  connector,  n-type,  p-type,  and  inter- 
connector  (exposed  to  ambient  as  shown  in  Fig.  1,  marked  with  (i)) 
are  subjected  to  convection  and  radiation  heat  transfer  boundary 
conditions. 

In  the  thermoelectric  material,  the  current  flow  and  energy 
transport  are  governed  by  the  continuity  of  current  density  and  the 
energy  equations  subjected  to  following  assumptions:  the 


materials  are  homogeneous  and  isotropic  with  temperature- 
dependent  properties;  the  thermoelectric  properties  such  as  elec¬ 
trical  resistivity  p,  Seebeck  coefficient  a ,  and  thermal  conductivity  k 
are  specified  in  polynomial  functions  of  temperature  (as  given  in 
Table  1  [15]);  the  thermal  and  electrical  contact  resistances  at  the 
interface  of  the  thermoelectric  and  copper  materials  are  neglected. 

The  system  of  partial  differential  equations  governing  the 
thermoelectric  phenomenon  in  thermoelectric  materials  under 
steady-state  conditions  for  the  current  flow  continuity  and  heat 
transport  [19]  are  written  as: 


V  J  =  0 


a) 


v(/<vr)  +  p]2 


Ty 


(V(x)j  + 


=  0. 


(2) 


In  Eq.  (2),  the  second,  third  and  fourth  terms  on  the  left-hand 
side  of  the  second  equation  represent  the  Joule  heating,  Peltier 
and  Thomson  effects,  respectively.  Although  the  Thomson  effect  is 
not  of  primary  importance  as  compared  to  Joule  and  Peltier  effects 


Table  1 

Grid  independence  study  for  composite  and  conventional  thermoelectric  devices  (grid  size  in  bold  face  is  chosen  for  further  simulations). 


Cells 

PoW 

%  of  error  |P°ne”  P°oldl  x  100 

I  Onew 

Qh  W 

%  of  error  IQhne™  <2holdl  x  100 
t^inew 

7]  % 

%  of  error  lVnew  ,,oldl  x  100 

Vnew 

Composite  TED  [at  Th  = 

550  K,  Tc  =  Tatm  =  Tsur  =  300  I<,  d  =  1  mm,  H=  10  mm,  h 

=  20  W/m2  K  and  RL  =  1.0  x  10"3  Q] 

11,6928 

0.6856 

13.2015 

5.1930 

208,568 

0.6843 

0.1802 

13.2492 

0.3604 

5.1650 

0.5426 

444,800 

0.6831 

0.1823 

13.2357 

0.1022 

5.1608 

0.0800 

877,500 

0.6815 

0.2235 

13.2105 

0.1908 

5.1591 

0.0326 

Conventional  TED  [at  Th 

=  550  K,  Tc  =  Tatm  =  Tsur  =  300  K,  H  = 

10  mm,  h  =  20  W/m2  K  and  RL  =  1.15  x  10"2  Q] 

102,076 

0.2007 

3.4341 

5.8443 

204,732 

0.2008 

0.0486 

3.4371 

0.0890 

5.8420 

0.5426 

393,600 

0.2007 

0.0429 

3.4360 

0.0318 

5.8413 

0.0800 

528,640 

0.2007 

0.0082 

3.4362 

0.0033 

5.8407 

0.0326 
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Table  2 

Comparison  of  thermoelectric  generator  parameters  Th  =  427  °C  and  Tc  =  27  °C. 


Quantity 

Analytical 

[21] 

ANSYS  (a) 
[22] 

ANSYS  (b) 
[22] 

Present 

(a) 

Present 

(b) 

Qh,  w 

13.04 

13.03 

11.07 

13.0 

10.72 

Po,W 

1.44 

1.43 

1.05 

1.43 

1.07 

ViPolQh),  % 

11 

11 

9.5 

11 

9.98 

J,  A 

19.2 

19.1 

16.4 

19.09 

16.5 

in  thermoelectric  devices,  it  has  been  considered  for  completeness 
of  analysis.  The  governing  equation  for  heat  transport  in  the 
connectors  and  inter-connectors  can  be  obtained  from  Eq.  (2)  by 
merely  removing  the  Peltier  and  Thomson  effects. 

Using  the  non-Ohmic  current- voltage  [20]  relationship,  the 
total  electric  potential  is  calculated  as: 


W  =  W0  +  Ws  =  -pj  -  aVT. 


(3) 


y0 

A^(RX  +tfL) 


,  dT 

and  -=■  =  0, 

a? 


(4) 


where  Vo  is  the  Seebeck  voltage.  The  Seebeck  voltage  is  the  open 
circuit  voltage  produced  when  a  temperature  differential  is 
imposed  across  the  junctions  of  the  device  at  a  no  load  condition. 
Furthermore,  V0  is  the  summation  of  the  Seebeck  potentials  at  the 
interfaces  of  the  semiconductor  and  conductor  materials  and  is 
expressed  as 


The  above  Seebeck  potentials,  with  reference  to  the  hot  side 
temperature,  are  evaluated  at  the  cold  side  of  the  interface  of  the 
semiconductor  and  inter-connector  materials. 

In  Eq.  (4),  R[  and  RL  are  the  total  internal  and  external  load 
resistances,  respectively.  The  R\  of  the  CTED  is  evaluated  as 


The  first  term  on  the  right-hand  side  of  Eq.  (3)  is  the  electrostatic 
potential  distribution  due  to  produced  electric  current  and  the 
second  term  represents  the  electric  potential  distribution  from  the 
Seebeck  effect  via  the  temperature  gradient  within  the  thermo¬ 
electric  material.  Since  the  temperature  distribution  and  current 
density  are  three-dimensional  (3D),  one  would  expect  the  voltage 
distribution  to  be  3D  as  well. 

The  associated  thermal  and  electrical  boundary  conditions  for  Eqs. 
( 1 )— (3 )  with  respect  to  the  geometry  as  shown  in  Fig.  1  are  as  follows. 

At  the  ‘in’  terminal: 


Current  density  n-type  leg  J  (A/cnT) 
-60  -45  -30  -15 


i  =  n,p,c,ic  1 


V,  J  ^ 


(6) 


where  V  is  the  volume  of  the  material. 
At  the  ‘out’  terminal: 


Fig.  2.  Comparison  of  efficiencies  for  example  1  in  Ref.  [19]  (a)  efficiency  of  n-type  and 
p-type  legs  (b)  module  efficiency  of  p-type  leg. 


Fig.  3.  Effects  of  load  resistance  on  (a)  electrical  current  and  (b)  Ohmic  voltage  for 
various  hot  wall  temperature  Th  values  [at  d  =  1  mm,  h  =  20  W  m-2  K1, 
TC  =  T(X  =  Tsur  =  300  K]. 
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V  =  0  and 


At  the  top  wall  of  the  upper  connector: 


where  ?  denotes  the  direction  normal  to  the  corresponding  wall  of 
(7)  the  element. 

The  power  output  from  the  TED  for  a  given  load  resistance  Ri 
and  heat  input  at  the  hot  wall  are  evaluated  as 


T  =  Th  and 


(8) 


Po  =  I2Rl,  Qh 


dT, 

Kc-^dAs 


(12) 


At  the  bottom  wall  of  the  lower  connector: 


and  the  thermal  conversion  efficiency  is  calculated  as 


T  =  Tc  and 


At  all  other  surfaces  exposed  to  surroundings: 


(9) 


V 


Po 

Qh' 


(13) 


q"  =  h(T  —  Too)  +  £<j(t4  -  7^ur). 


(10) 


Here,  in  Eq.  (10),  the  first  and  the  second  terms  on  the  right  hand 
side  represent  the  effects  of  convention  heat  transfer  and  radiation 
on  TED,  respectively.  The  emissivities  t  for  the  semiconductor  and 
the  highly  polished  copper  conductor  are  taken  as  0.45  and  0.03, 
respectively.  The  ambient  fluid  (Too)  and  the  surrounding  (Tsu r) 
temperatures  are  kept  constant. 

At  the  interface  between  the  semiconductor  and  conductor 
materials,  the  continuity  of  temperature,  current  density  and  the 
heat  flux  conditions  are  imposed  and  are  written  as 


^c,ic  —  ^n,p?  Jc,ic  —  Jn,p:>  cUld 


dTcM 

d? 


^n,p  dTn,p 

Kc,ic  d? 


(ID 


3.  Numerical  solution  procedure  and  validation 

The  numerical  simulations  are  performed  in  a  FLUENT-UDS 
(User  Defined  Scalar)  environment  using  the  finite  volume 
formulation  of  the  governing  Eqs.  (1)  and  (2),  and  the  constitutive 
relation  (Eq.  (3))  along  with  the  associated  boundary  conditions 
(Eqs.  (4)  and  (7)).  FLUENT  has  well-implemented  fluid  and  heat 
transfer  models,  allowing  us  to  numerically  investigate  the  thermo¬ 
fluid-electrical  coupling  mechanisms  in  TEDs.  Additionally, 
numerical  solutions  for  a  thermoelectric  generator  have  been 
demonstrated  using  the  UDS  package  in  FLUENT  [13]. 

The  geometric  models  and  mesh  are  generated  in  Gambit  2.4. 
The  Ohmic  heating,  Peltier  and  Thomson  effects  are  modeled  as 
source  terms  in  the  energy  equation  (Eq.  (2)).  The  Ohmic  and 


Fig.  4.  Effects  of  load  resistance  on  (a)  electrical  power  output  (b)  heat  input  and  (c)  thermal  efficiency  for  various  hot  wall  temperature  Th  values  [at  d  =  1  mm,  h  =  20  W  m  2  K  \ 
Tc  =  Too  =  TSUr  -  300  K]. 
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Seebeck  potential  distributions  (Eq.  (3))  are  calculated  using  UDS 
fields.  The  spatial  discretization  of  the  diffusion  term  is  done  with 
a  power  law  scheme.  For  a  given  load  resistance  Ri,  the  electric 
current  (as  expressed  in  (Eq.  (4)))  is  evaluated  based  on  the  open- 
circuit  Seebeck  voltage  (Eq.  (5))  generated  via  the  temperature 
differential  applied  at  hot  and  cold  junctions. 

The  convergence  criteria  (the  absolute  residual  difference 
between  consecutive  iterations  over  a  domain)  for  Ohmic  and 
Seebeck  potentials,  and  energy  are  set  as  10  10  and  1CT12, 
respectively.  Grid  independence  tests  are  performed  for  both 
composite  and  conventional  TEDs  (as  shown  in  Fig.  1)  and  are 
depicted  in  Table  1.  Grid  sizes  with  208,568  and  204,732  cells  are 


chosen  for  the  analysis  of  composite  and  conventional  TEDs, 
respectively.  The  orthogonal  and  non-uniform  grid  is  employed  for 
succinctness  and  the  computational  mesh  used  for  simulations  is 
shown  in  Fig.  1. 

The  three-dimensional  thermoelectric  model  has  been  validated 
with  published  results  of  a  thermoelectric  generator  for  the  cases  of 
(a)  constant  properties  [21]  and  (b)  temperature-dependent 
properties  [22]  as  given  in  Table  2.  Fig.  2  also  depicts  the  valida¬ 
tion  of  the  present  numerical  code  with  one  of  the  examples  from 
[19]  in  terms  of  individual  leg  and  module  efficiencies.  The  results 
are  in  good  agreement  with  published  analytical  and  numerical 
results. 
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Composite  TED  Conventional  TED 

Fig.  5.  The  contours  of  (a)  temperature  (b)  Seebeck  and  (c)  Ohmic  voltages  at  optimum  load  resistance  values  for  both  composite  and  conventional  TEDs  [at  Th  =  550  K,  d  =  1  mm, 
h  =  20  W  m"2  K1,  TC  =  T„  =  Tsur  -  300  K]. 
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4.  Results  and  discussion 

The  numerical  simulations  are  conducted  on  composite  ther¬ 
moelectric  devices  (CTEDs)  as  well  as  a  convectional  TED  made 
with  n-type  (75%  Bi2Te  25%  Bi2Se),  p-type  (25%  Bi2Te3  75%  Sb2Te3 
(1.75%  excess  Se)),  and  copper  materials.  The  connector  has 
a  geometry  of  length  by  depth  by  height  of  35  x  5  x  1  mm3.  The 
height  of  the  thermoelectric  leg  is  H  =  10  mm  and  the  cross- 
sectional  area  of  the  n-  or  p-type  leg  is  5  x  5  mm2.  The 
distance  between  the  legs  is  I2  =  5  mm.  The  cold  wall,  ambient 
and  surrounding  temperatures  (Tc  =  Too  =  Tsur  =  300  I<),  and  the 
geometric  dimensions  as  given  above  are  invariant  in  the  present 
study.  The  influence  of  the  hot  wall  temperature  350  <  Th 
(I<)  <  550,  the  load  resistance  10~4  <  Pl  (Q)  <  10-1,  the  semi¬ 
conductor  thickness  0.25  <  d  (mm)  <  4  and  the  convection  heat 
transfer  coefficient  0  <  h  (W  m_1  K_1)  <  103  on  the  thermo¬ 
electric  performance  of  a  CTED  compared  with  geometrically 
equivalent  conventional  TED’s  values  has  been  investigated  in 
detail. 

Considering  a  geometrically  equivalent  CTED  with  d  =  1  mm 
and  a  conventional  TED  for  a  given  h  =  20  W  m-1  K-1  value,  Fig.  3 
shows  the  variation  of  the  electric  current  I  and  the  Ohmic  voltage 
drop  Vo  while  Fig.  4  depicts  the  heat  input  Qh,  power  output  P0  and 
the  thermal  conversion  efficiency  77  with  respect  to  load  resistance 
Ri  for  different  values  of  hot  wall  temperatures  Th.  In  these  figures, 
the  solid  lines  with  filled  circles  represent  the  results  for  a  CTED  and 
the  dotted  lines  with  empty  circles  describe  conventional  TED 
values.  From  Figs.  3  and  4,  it  is  observed  that  for  a  given  PL  value,  an 
increase  in  Th  results  in  an  increase  in  /,  V0,  P0,  Qh  and  77  values. 
However,  at  a  given  Th,  an  increase  in  Ri  decreases  both  I  and  Vo 
values,  with  these  values  eventually  approaching  zero  at  larger  PL 
value.  This  indicates  the  TED  produces  zero  Po  at  large  Ri  value  even 
though  the  device  experiences  a  temperature  differential  Th  -  Tc. 
The  similar  behavior  is  also  observed  in  the  conventional  TED.  For 
a  given  PL,  the  CTED’s  I  value  is  higher  as  compared  to  the 
conventional  TED.  On  the  other  hand,  the  opposite  behavior  is 
observed  for  Vo.  This  is  due  to  the  total  internal  resistance  R[ 
(summation  of  copper  and  semiconductor  materials,  i.e. 
Pc  +  Pn  +  Pp)  of  the  CTED  being  less  as  compared  to  the  conven¬ 
tional  TED.  It  is  noticed  that  the  difference  between  the  ordinate 
values  of  the  CTED  and  the  conventional  TED  decreases  with  PL  in 
case  of  /,  whereas  the  similar  behavior  has  not  been  seen  for  Vo. 
Moreover,  when  Ri  <  2  x  10  2,  there  is  a  substantial  change  in  I  and 
Vo  values. 

In  Fig.  4a  and  c,  for  a  given  Th,  it  is  observed  that  there  exists 
a  maximum  Po  when  Ri  =  R[  and  an  optimum  77  at  a  particular 
resistance  PL  known  as  the  optimum  load  resistance  PoptmL-  For 
the  CTED,  the  R[  values  (in  Q,  as  calculated  from  Eq.  (6))  evaluated 
at  rh  =  350,  450  and  550  K  are  1.855  x  10“3,  2.217  x  10“3,  and 
2.254  x  10-3,  respectively;  for  the  conventional  TED,  the  values 
are  8.645  x  10“3,  9.991  x  10“3,  and  1.1135  x  10“2.  The  RoptmL 
increases  with  an  increase  in  Th  values  while  all  other  parameters 
are  kept  constant.  The  PoptmL  values  for  the  CTED  are  one  order  of 
magnitude  less  than  the  conventional  TED,  and  interestingly 
enough,  the  CTED’s  behavior  is  extremely  similar  to  that  of  the 
conventional  TED.  However,  greater  Po  values  and  lesser  PoptmL 
values  are  predicted  for  the  CTED  when  compared  to  the 
conventional  TED.  It  is  also  observed  that  the  optimum  77  values  of 
the  CTED  are  24.8%,  26.2%  and  29.9%  higher  than  the  conventional 
TED  at  Th  =  550  K,  450  K  and  350  K,  respectively.  This  trend  is  due 
to  Po  and  Qh  being  functions  of  Pi;  as  Pi  decreases  for  the  same  Th 
value,  Qh  increases,  and  the  inverse  is  true;  Qh  decreases  with  an 
increase  in  Pl  value.  Moreover,  within  the  specified  Pl  range,  the 
CTED  is  subjected  to  higher  Th  values  and  shows  a  dramatic 
change  in  Qh  predictions.  As  shown  Fig.  4b,  the  CTED  exhibits 


nearly  a  three-fold  increase  in  Qh  values  when  compared  to  the 
conventional  TED. 

For  a  given  Th  =  550  K  and  h  =  20  W  m-1  K_1,  the  three- 
dimensional  (3D)  distributions  of  temperatures,  Seebeck  and 
Ohmic  potentials  for  both  composite  (d  =  1  mm)  and  conven¬ 
tional  TEDs  are  plotted  in  Fig.  5.  Here,  the  PL  is  taken  at  the 
optimum  77  value  (Fig.  4)  for  the  corresponding  TED  case.  For 
a  given  Th  -  Tc,  irrespective  of  the  TED  case,  the  temperature  has 
to  flow  from  the  hot  side  surface  to  the  cold  side  surface.  It  is 
observed  that  the  temperature  and  voltage  gradients  in  the 
semiconductor  material  of  the  CTED  are  higher  compared  to  the 
conventional  TED.  It  is  also  observed  that  the  temperature  and 
Ohmic  voltage  gradients  are  negligible  in  the  conductor  material 
as  compared  to  the  semiconductor  material.  This  is  attributed  to 
the  high  thermal  and  electrical  conductivity  of  the  conductor 
material.  It  can  be  seen  that  the  presence  of  lower  a  convection 
heat  transfer  coefficient  results  in  minimal  nonconformity  in  the 
temperature  and  voltage  distributions  of  leg.  Due  to  larger  Pj  of 
the  conventional  TED,  the  Ohmic  voltage  drop  and  the  total  built- 
in  Seebeck  voltage  values  are  higher  when  compared  to  the  CTED 
values. 


Fig.  6.  The  variation  of  (a)  electrical  current  and  (b)  Ohmic  voltage  with  load  resistance 
for  different  semiconductor  thickness  d  values  [at  at  Th  =  450  K,  h  =  20  W  m-2  K-1, 
TC  =  T„  =  Tsur  =  300  K]. 
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Figs.  6  and  7  show  the  electrical  and  thermal  characteristics  (/, 
Vo,  Po,  Qh,  and  7j)  of  a  CTED  for  different  semiconductor  thickness 
d  and  load  resistance  Ri  values  with  conditions  of  Th  =  450  K  and 
h  =  20  W  m2  K-1.  The  CTED  qualitatively  shows  similar  behavior 
for  a  given  d  value  as  illustrated  in  Figs.  4  and  5.  Furthermore,  for 
a  specified  PL  value,  an  increase  in  d  results  in  a  decrease  in  I  and 
an  increase  in  Vo-  There  is  a  significant  change  in  I  when 
Ri  <  2  x  10-3;  however,  a  moderate  change  in  Vo  for  the  range 
10-4  <  Ri  <  10-1  is  observed.  This  behavior  is  explained  by  the 
cumulative  effect  of  R[  and  Ri  on  I  for  a  given  Th  -  Tc  value.  For 
instance,  the  contours  of  the  Seebeck  Vs  and  Ohmic  V0  potentials 
for  various  d  values  at  P0ptmL  are  shown  in  Fig.  8.  The  percentage 
change  in  the  total  absolute  value  of  Vs  (as  depicted  in  Fig.  8)  is 
minimal  with  d  when  compared  to  the  change  in  Vo  values.  This  is 
due  to  the  temperature-dependence  of  the  Seebeck  coefficient 
and  the  minimal  temperature  drop  across  the  inter-connector 
material. 

From  Fig.  7c,  it  is  seen  that  as  d  increases  the  optimum  tj  value 
increases,  and  it  reaches  a  maximum  value  at  d  =  1  mm.  After  that, 
the  optimum  rj  starts  decreasing  with  an  increase  in  d.  Interestingly, 
the  optimum  rj  value  of  the  conventional  TED  is  higher  than  CTED 
case  when  d  =  4  mm.  The  interplay  of  geometric  dimensions,  heat 
loss  to  surroundings  and  the  boundary  conditions  at  the  hot  and 
cold  surfaces  explains  the  effect  on  the  optimum  7]  values.  From 
Fig.  7a  and  c,  both  Po  and  Qh  substantially  increase  with  a  decrease 
in  d  values.  Additionally,  the  P0ptmL  value  increases  with  an  increase 
in  d.  Irrespective  of  the  d  value,  beyond  Ri  =  10-1,  the  Po  reaches 
a  single  value  which  represents  the  conventional  TED  value  (as 
shown  in  Fig.  7a).  The  Qh  as  seen  in  Fig.  7b  shows  a  higher  value  at 


d  =  0.25  mm  for  a  given  Ri  value  as  compared  to  other  d  values. 
Moreover,  as  semiconductor  thickness  d  doubles,  heat  input  Qh 
approximately  reduces  to  half  value  and  vice-versa.  This  increase  in 
Qh  with  reduction  in  d  is  due  to  the  rise  in  the  effective  thermal 
conductance  of  the  leg  while  maintaining  a  fixed  Th  -  Tc  value.  The 
temperature  distribution  in  the  CTED  for  different  d  values  is 
depicted  in  Fig.  8.  It  is  noticed  that  the  temperature  gradients  in  the 
semiconductor  increases  with  a  decrease  in  d  and  vice-versa.  From 
Fig.  8,  it  is  also  observed  that  at  a  given  h  value,  as  d  increases,  the 
distributions  of  temperature  and  voltages  become  more  three- 
dimensional. 

For  different  values  of  d  and  Th,  the  effects  of  the  convection  heat 
transfer  coefficient  h  (represents  the  heat  loss  from  the  surfaces 
exposed  to  the  surroundings)  on  the  thermoelectric  characteristics 
of  the  CTED  are  shown  in  Figs.  9  and  10.  Here,  in  the  simulations, 
the  Ri  is  set  equal  to  R[  (as  depicted  in  Fig.  9c).  For  given  d  and  Th 
values,  the  variations  in  I,  Vo,  and  Po  are  negligible  with  an  increase 
in  h  from  5  to  100;  however,  when  h  >  100,  a  noticeable  change  is 
observed  in  these  values.  Contrary  to  this  trend,  Fig.  10b  and  c 
illustrate  a  significant  change  in  Qh  and  rj  values  for  the  same  range 
of  h  values.  In  order  to  maintain  the  constant  Th  -  Tc  condition,  the 
Qh  value  is  substantially  increased,  and  hence  rj  is  reduced  with  an 
increase  in  h  value.  It  is  also  observed  that  the  influence  of  h  on  R[  is 
strong  at  higher  Th  values  as  shown  in  Fig.  9c.  For  brevity,  when 
h  =  0  (no  heat  loss)  condition,  the  performance  parameter  (I,  Vo,  R[, 
P0,  Qh,  and  77)  values  are  given  in  these  figures  for  the  corresponding 
d  and  Th  values. 

To  put  these  findings  into  perspective  of  practical  application, 
a  simple  cost  analysis  is  presented.  The  percent  reduction  of 


Fig.  7.  The  variation  of  (a)  electrical  power  output  (b)  heat  input  and  (c)  thermal  efficiency  with  load  resistance  for  different  semiconductor  thickness  d  values  [at  Th  =  450  K, 
h  =  20  W  m"2  K1,  TC  =  T„  =  Tsur  -  300  K]. 
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Fig.  8.  The  distributions  of  temperature,  Seebeck  and  Ohmic  voltages  for  different  semiconductor  thickness  d  (mm)  (a)  0.25  (b)  1  (c)  4  and  (d)  5  (conventional  TED)  at  optimum  load 
resistance  values  [at  Th  =  450  K,  h  =  20  W  m"2  K1,  Tc  =  Tx  =  Tsur  =  300  K]. 


material  usage  and  associated  costs  of  a  composite  leg  are 
compared  to  a  conventional  leg  as  shown  in  Table  3.  Based  upon 
laboratory  experiments,  the  price  of  commercially  available 
bismuth-telluride  semiconducting  material  is  approximately 
$0.72  per  mm3.  The  corresponding  high  thermal  conductivity 
oxygen  free  copper  materials  cost  $0.0069  per  mm3.  For  instance, 
the  cost  of  a  composite  leg  (Fig.  la)  with  d  =  0.5  mm,  L  =  10  mm 
and  a  cross-section  of  5  mm  x  5  mm,  (includes  thermoelectric 
material  and  inter-connector)  would  be  $19.55,  as  compared  to  the 
conventional  element  with  a  cost  of  $180.  This  reflects  a  rare-earth 
element  material  usage  reduction  of  90%  and  a  cost  reduction  of 


89.14%.  Although  the  cost  and  rare-earth  element  material  usage  is 
reduced,  this  data  does  not  reflect  the  extra  processing  needed  to 
construct  the  composite  leg:  polishing  of  the  copper  surfaces 
which  will  be  in  contact  with  the  bismuth-telluride;  the  nickel 
coating  of  the  copper  to  prevent  oxidation  and  rapid  copper 
diffusion  into  the  bismuth-telluride;  the  diffusion  barrier  depo¬ 
sition  on  top  of  the  nickel  coating  to  ensure  prevention  of  nickel 
and  copper  diffusion  into  the  semiconductor;  the  diffusion 
bonding  process  itself.  All  of  these  additional  manufacturing  steps 
do  not  have  a  quantifiable  cost  associated  with  the  production  of 
the  composite  leg. 


a  100 


Fig.  9.  The  influence  of  convection  heat  transfer  coefficient  on  (a)  electrical  current  (b)  Ohmic  voltage  and  (c)  internal  resistance  for  various  hot  wall  temperature  Th  and  semi¬ 
conductor  thickness  d  values  [at  RL  =  Rit  Tc  =  =  Tsur  =  300  K]. 
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Fig.  10.  The  influence  of  convection  heat  transfer  coefficient  h  on  (a)  electrical  power  output  (b)  heat  input  and  (c)  thermal  efficiency  for  various  hot  wall  temperature  Th  and 
semiconductor  thickness  d  values  [at  RL  =  R[,  TC  =  TX  =  Tsur  =  300  K]. 
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Table  3 

Quantitative  comparison  of  composite  leg  material  cost  and  usage  with  conven¬ 
tional  leg  (H  =  10  mm,  W=  5  mm  and  D  =  5  mm). 


Semiconductor 
element  size  (d,  mm) 

Cost  of  composite 
leg  (in  US  dollars) 

%  of  reduction 
in  material  cost 

%  of  reduction  in 
material  usage 

0.25 

10.64 

94.09 

95 

0.5 

19.55 

89.14 

90 

1 

37.38 

79.23 

80 

2 

73.04 

59.43 

60 

4 

144.35 

19.81 

20 

5.  Conclusions 

Using  the  FLUENT-UDS  environment,  the  thermoelectric 
coupling  mechanisms  of  Joule  heating,  Peltier  cooling,  and  Thom¬ 
son  cooling/heating  in  a  thermoelectric  material  with  temperature- 
dependent  properties  have  been  implemented.  Oriented  toward 
improving  a  thermoelectric  device’s  (TED)  power  output  and  also 
minimizing  the  usage  of  rare-earth  materials  (n  or  p-type)  in  waste 
heat  recovery  applications,  a  novel  composite  thermoelectric 
device  (CTED)  has  been  proposed  without  compromising  its 
thermal  conversion  efficiency. 

The  thermoelectric  performance  of  a  CTED  compared  with 
a  conventional  TED  in  terms  of  electric  current  /,  Ohmic  Vo  and 
Seebeck  Vs  voltages,  power  output  P0,  heat  input  Qh  and  thermal 
conversion  efficiency  r/  for  different  values  of  hot  surface  temper¬ 
ature  Th,  load  resistance  Ri,  semiconductor  thickness  d,  and 
convection  heat  transfer  coefficient  h  are  investigated  in  detailed 
using  numerical  simulations. 

The  difference  in  power  outputs  of  composite  and  conventional 
TEDs  is  decreased  with  an  increase  in  Ri  while  all  other  parameters 
are  kept  constant.  It  is  demonstrated  that  P0  reaches  a  maximum  at 
Ri  equaling  the  total  internal  resistance  Pi  and  tj  attains  an  optimum 
at  an  optimum  load  resistance  PoptmL.  For  the  CTED,  the  PoptmL  is 
lower  than  the  conventional  TED  value.  For  both  TEDs,  a  marginal 
increase  in  PoptmL  value  with  an  increment  in  Th  is  noticed.  It  is  also 
observed  that  a  CTED  with  d  =  1  mm  shows  a  24.8%,  26.2%,  and 
29.9%  increase  in  ri  as  compared  to  a  geometrically  equivalent 
conventional  TED  when  Th  =  550  K,  450  K  and  350  K,  respectively. 
At  the  minimum  stated  temperature,  the  CTED  showed  a  near 
three-fold  increase  in  Qh- 

For  the  CTED,  it  is  noticed  that  as  d  increases  the  optimum  7] 
value  increases,  and  it  reaches  a  maximum  value  at  d  =  1  mm.  After 
that,  the  optimum  r/  starts  declining  with  an  increase  in  d.  On  the 
other  hand,  both  Po  and  Qh  substantially  increase  with  a  reduction 
in  the  d  value.  At  higher  Ri  values,  irrespective  of  the  d  value,  the  Po 
value  reaches  that  of  a  conventional  TED.  Moreover,  as  d  is  reduced 
by  half,  the  Qh  approximately  doubles. 

For  a  given  CTED  configuration  d  and  temperature  differential 
Fh  -  Tc,  the  convection  heat  transfer  coefficient  h  has  a  negligible 
effect  on  Po.  Elowever,  h  has  a  marginal  influence  on  both  Qh  and  rj 
values. 

Finally,  it  is  proposed  that  building  state-of-the-art  novel  CTEDs 
with  a  minimum  amount  of  rare-earth  metals  reduces  cost  without 
compromising  thermal  conversion  efficiencies  and  performance, 


greatly  helping  with  the  production  of  clean  energy  from  waste 
heat  resources. 
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